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Abstract 

We compute the Casimir Energy of a spherical region using a Surface Impedance 
approach. We characterize the Surface Impedance of the boundary using plasma model. 
Exact analytical formulae are obtained by means of the zeta function regularization 
method and all the divergencies are explicitly computed. We find that it is possible 
to have negative Casimir energy for some range of values of the relevant parameter 
y a = uj v aJT[i. Limits of applicability of the model are discussed. 
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1 Introduction 

The interest in the Casimir energy in spherical bodies is enormous. The first attempt of 
Casimir to furnish a model for the electron pQ in which an attractive Casimir force would 
balance the electrostatic self-repulsion showed to be unsuccesful after the finding of Boyer 
[2] that the Casimir force in a sphere is repulsive. The interest in this configuration is 
still very strong and the use of advanced computational tools has simplified very much the 
calculation (see [3] and references therein). Here we make another small step, aimed at 
further simplifying the procedure, implementing the Surface Impedance approach in this 
kind of problems. 

Surface Impedance (SI) can result very useful in modeling the behavior of a non ideal 
surface. It is defined through the equation [H [61 [7] : 

B tan \ a = Z(nxB)\ a (1) 
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with n the outward normal to the surface. This formula relates the tangential fields outside 
the material surface and all the characteristics of the material are taken into account 
through the values of Z. Eq. (pQ) can be seen as the defining equation for Z and it can be 
applied to arbitrary materials even when a description in terms of dielectric permittivity 
cannot be given 0[9]. In particular it can result very useful in describing the transition 
from attractive to repulsive behavior |10]lll| when computing Casimir energy [121 113] and, 
in this respect, this paper can be seen as the natural prosecution of [llj. (As far as we 
know the first time SI has been used in connection with Casimir energy was in [7] ). 

The approach is very general and powerful and, in our opinion, it can treat, in a 
relatively easy manner, problems with spherical boundaries: a dielectric sphere |14j . the 
bag model of quantum chromodynamics [151 [T6] , quantum gravity and cosmology [17] [T8] , 
boundaries induced by topological parameters [19] etc. Obviously the possibility of using 
an equation like (pQ) must be treated carefully case by case, here we make a first attempt 
having in mind "standard cases" like electromagnetic field in a sphere and the MIT bag 
model of QCD [16] 120] . It is important to stress from the very beginning that within 
this model the renormalization procedure is not completely settled and the model is still 
physically unsatisfactory [211 [3]. To regularize the Casimir energy we will use the zeta 
function regularization p2| [23l [2^1 [25] [26| [3] . 

The paper is organized as follows: in Sec. II the zeta function regularization approach 
is shortly revised for the convenience of the readers. In Sec. Ill the relevant formulae are 
derived and the Casimir energy is computed, finally in Sec. IV we analyse the results. 



2 Casimir Energy and Zeta function regularization 

The Casimir energy is the vacuum expectation value of the Hamiltonian operator on the 
ground state: 

E Cas =< 0\H\0 >=i^(^j) (2) 

J 

where Ej are the energy eigenvalues labelled by some general index J. The sum is, in 
general, divergent and regularization is necessary. In the £ function regularization scheme 
one defines a new quantity (h(s) as 

E Cas = Mm J^-J2( E jV S =-- ^y s+1 (H(s) (3) 

s— >— 1/2 z, s— 1/2 

where Ch( s ) is the £ function relative to the differential operator connected to the operator 
H. The parameter /i is an arbitrary parameter introduced for dimensional reasons, it will 
disappear on removing the regularization in the limit s — > —1/2. In general no explicit 
expression for the eigenvalues Ej exists, however one can use the argument principle to 
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represent the sum over the eigenvalues as a contour integral in the complex plane \27\ [28] : 

9{a k ) ~ 9(bm) = ^-j> 9{z)dlog(A{z)) (4) 

km J 1 

where 7 is a closed contour containing all the zeros a n and poles b n of the function A(z) 
assumed analytic in 7 and g{z) is some analytic function inside 7. Thus if A(z) is such 
that A(u n ) = 0, u n being the eigenvalues of our problem, and has no pole, the sum over 
uj n can be obtained as a contour integral. In this case A{z) is called, for obvious reasons, 
the mode-generating function. In conclusion can be written as a contour integral of 
some mode-generating function [29 : 

(H = Y^-^g(z)dlog(Aj(z)) (5) 
where the sum over J takes into account possible degeneracy of the eigenvalues. 



3 The Casimir Energy for a Sphere 

In the following we will concentrate on the case of an electromagnetic field in a sphere of 
radius a. We will characterize the boundary by means of its Surface Impedence. Assuming 
a time dependence e~ lult the electric and magnetic fields in the interior of the sphere can 
be written in the form [6j O [301 [31] 



E = E ik aTE [™jAkr)Y lm (6, <)>) + (krj u (kr)Yn x X lm ] + a™j u (kr)X, 



1=1 

00 



1 in 



k 



H = Y,— \ aTE Mkr)X 



^lm ~ 1- 



rk 



[ihj u (kr)Y lm (9,(f>) + {krj u {kr))'ii x X, 



l.rn 



where Y/ m and X; m are the scalar and the vectorial spherical harmonics respectively, k = 
yfejhjj, ji/(x) = \/^Ji+i/2(x) are the spherical Bessel functions [H [33], and (xf(x))' = 

Imposing boundary conditions: Eq. dl|), we find the equations for the TE and TM 
modes: 



A™(x) 



— (kaj u (ka))' - Zj u (ka) 

-Z^- {kaj u {ka))' + j u {ka) 
ka 



a TE = 



a™ = 0, 



(6) 
(7) 
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where and A^ M are our mode-generating functions. In this case can be given in 
the following form \2b\ [3] : 



n=0 1=0 ^ ' 



m~ 



(8) 



where are the eigenmodes, H = 1, c = 1 is assumed. We introduced a (fictitious) mass 
parameter that we will let go to zero in the end of the calculation so to avoid some problems 
with the integral representation of the C function. Using formula ([5]), with a suitable choice 
of the contour 7, (jy can be written [26] 



0) 



with 1/ = 1 + 1/2. Shifting the integration contour along the imaginary axis and using 
ji(x) = \p^.Jv{x) and J v (ix) = e mv J v {— ix) and J u (ix) = e w 2/ l/ (x), ^(x) being the 
modified Bessel functions [33], we obtain the following expression valid in the strip 1/2 < 
»(a) < 1. 



5> 

^=3/2 



sin(7rs) 

7T 

sin(7rs) 

u=3/2 

t TE (s) + (™(s) 



dk [k 2 - m 2 ] 



ma/v 



dy 



- s d_ 

dk 



m 



log[k- 2v A^(ka)Ai; M (ka) 



d_ 

dy L 



log(y- u A TE (yu))+log[y- u A 



where, after rotation to the imaginary axes, the mode-generating functions are given by: 



A 



TM, 



X) 



h(x) 



Z(ix/a) 
2x 



— — Z(ix/a)x 



- Z{ix / a)i v {x) 
+ xl v {x) 



(12) 
(13) 



It is clear, now, the role played by the mass m. Indeed representation (jlpp is defined for no 
values of s if m = 0. The procedure can be modified for m = but it is more difficult [34J. 
Unfortunately we need the zeta function computed to the left of the strip 1/2 < 3?(s) < 1. 
Now, the restriction 1/2 < 3?(s) is due to the behavior of integrand as k — > 00. The general 
technique [26j to overcome this problem is to add and subtract the asymptotic term 



log ( y- v K TE {yv) ) + log ( y-»A 1M (yu 



-uXTM, 



to the integrand so to move the strip of convergence to the left. If we are able to compute 
analytically the asymptotic term alone and at the same time to treat the remaining term 



at least numerically (but, in general, it can be made very small by considering higher and 
higher terms in the asymptotic expansion) we can obtain the required analytical continua- 
tion. 

To be more specific, let us define Cn(s) by the following equality 



(ciroo - cr m (*)) + Cu ym {s) =■■ cn(s) + cr»- 



(14) 
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In the following we will give a representation of C,% v (s) in terms of known functions 
and valid in the region of interest of the complex plane. Obviously we have to choose an 
expression for Z to characterize the properties of the medium. We will use the plasma 
model; thus 

: yv} = y = . y 
V s u + y 2 ' x *« 

with xg u = \J b 2 + y 2 , 5 U = and y a = a^fejlujp, uj p being the plasma frequency of the 
material constituting the surface. At this point the procedure is quite standard j22j tZS 
261 [3] we need the asymptotic values of A TE and A™ for v — > oo with kjv fixed. To this 
extent we use the uniform asymptotic expansions of I v (x) and I u (x) [32| [33] : 



with 



Iv(vx) 

i v (vx) 



yfhuj (1 



,2^1/4 



u k (t) 



fc=l 



y+y 



2\l/4 



V2 



k=l 



Vk(t) 



(15) 
(16) 



u (t) 
vo(t) 



1 



^, and rj = yjl + y 2 + In ; 

i 2 d_ t 2 ) i r* 

1, andu fc+1 (t) = -i— -<(*) + - / cfe(l-5z 2 H(z);/c = 0,l,2. 

2 o JO 



1, and Ufc(t) = ^^(t) — t(l — t 2 



Uk-i(t) + i«Jsj_i(t) 



;k = 0,1,2... 



3.1 TE-Modes 

Inserting the asymptotic expansions of Ij, and I u in ()10p we obtain: 



Msym 



sin(7rs) 



7T 



E 



i>=3/2 

-: Cq + Ci + C2 



ma/v 



yv 



m 



d_ 

dy 



log ( y-"k TE (yv 



(17) 
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with (in the following we perform our calculations up to n max = 4) 



log[y-»A TE (yv) 



log ^ y~ v 



vr) 



V2 



.2\l/4 



k=l 



so that 



sin(7rs) 



7T 



Ci 



c 2 



sin(7rs) 



7T 



■ E " 

^=3/2 

CO 

E " 

i/=3/2 



ma/v 



dy 



{-) - m 



d_ 

dy 



log 



VT] 



V2 



(1 + ^) 



2\V4 



^2/ 



ma/v 



yv\- 
a J 



m 



9 i r 1 i 



sinM g ^ d J^ 2 _ m 2]- S ( 



is=3/2 



a ,,a , E\lM j_ ju A 3(j/) , At(y) 

A)U/J H \- ,, 2 + — 3 1 -4— 



where the functions A^y), (i = 0,4) are given in the appendix A. 

For the sake of clarity the calculations are developed in Appendix B and here we report 
the results only: 



Co = 
Ci = 



48a 



X (2 S -2,|)sin(7r S )r( S -i)r(- S ) 



7T 



3/2 



sin(Trg) 2s 
O2 = a 



7T 



C ( 2s - 1, f) /™( a ) + C f 2a, ^ *(,) + C + 1, f) 



C^ + 2,|j/ 3 ^( S ) + c(2 S + 3,^/r( S ) 

with the ff E {s) for (i = 0,4) given in appendix B. 
expanding around s = — 1/2 we finally obtain: 



0.0353 



C£f m (-l/2) = 0.107 + 0.229y a + 0.127yi + 0.0706 log(a) + ^ I ^ + 



(18) 



2 / 0.0624 - 0.0398 log (a) 
' Ja V s + 1/2 



- 0.0398 log 2 (a) +0.125 log(a) - 



0.0199 ' 

(s + 1/2) 2 . 
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3.2 TM-Modes 

In the same manner we have 

c T T(s) = ^ E v d y 



TV 



i/=3/2 

=: A) + £>i+£>2 



majv 



m 



log (y-"A™(yi/)) 



(19) 



with 



logfy-'A^V) 



2\l/4 



fc=i 



. / 1 2 \ / nmax Nil 

H^Hs^HJl (20) 



and 



Di 



sin(7rs) 



7T 



i/=3/2 



majv 



yv 



* v % Jrna/u « > 



£>2 = 



Sin 7TS 



m 



m 



m 



9 1 



2\V4 



Co 



5- 1 



v=3/2 



Bi(y) , 5 2 (y) , ^(y) , B 4 (y) 



+ 



+ 



+ 



with the functions Bi, [i = 1,4) given in the appendix A. We find 
A) = Co 

a 2s C (2s - 1, |) (sin(^)r (± - s) Y{s) + vr 3 / 2 ) 



D, = — 



sin(7rs) 



,2s 



7T 



2vr 3 / 2 



C ( 28,$ ) "to + C ( 2* + 1,^ ] # M (s) + C ( 2s + 2, ^ ) 



TM 1 



cTE, 



si 2s + 3,^) // -»(.,) 



with the/?" M (s), (z = 1,4) defined in appendix B. Expanding aroud s = —1/2 we find: 



c 



asym 
TM 



-{-0.170 + 0.266y 2 - 0.317^ - 0.05541og(a) - 



0.0277 
s + l/2 



+ 



vl 



0.130 - 0.119 log(a) 



s + l/2 



0.119 log 2 (a) + 0.259 log(a) 



y 4 J 0.130 log(a) + 



0.065 
s + l/2 



)} 



0.0597 



+ 
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Thus the final result is (after having reintroduced the parameter fj,): 

E C as = -{-0.063 + 0.229y a + 0.393^ - 0.317^ + 0.015 log(//a) + °; ^ + 
a I s + l/2 

9 /0.192 - 0.159 log(ua) , . „ , , 0.0796 \ 

Vl ( s + 1/2 m ~ 0.159 log^a) + 0.3841ogM - j + 

^(o.l301og(H + ^^)} (21) 



4 Discussion 

Formula (|21 j) has many interesting features. First the structure of the divergencies: to- 
gether with the usual terms depending on s+ */ 2 an d logC/- 4 * 1 ) we nn d a new one: a second 
order pole in (s + l/2) and the log 2 (//a) term. This terms cannot be eliminated com- 
puting the Principal Part of the zeta functions as usually done within the zeta function 
regularization |35j . What we expect, in general, is that these divergencies cancel when 
exterior field modes are included [361 [3] . However the peculiar (and interesting) feature of 
the SI approach is exactly the fact that one can perform all the calculation without any 
reference to exterior modes. So it would be desirable that surface terms renormalize the 
occurring divergencies but, on the other side, it cannot be excluded that the renormal- 
ization procedure will induce further finite contributions to the calculated ones. In this 
respect the situation is similar to that of a shell of finite thickness in which one cannot 
give a satisfactory interpretation of the vacuum energy [3], and further investigations are 
necessary. 

More interesting, from the physical point of view, results the finite part of the energy: 

Ecas = ^{-°- 063 + °- 229 2/a + 0-393y 2 - 0.317^} (22) 

Indeed it turns into negative values both for small, y a < 0.2, and large, y a > 1.3, values of 
y a . In particular it is negative for y a = 0, but, it corresponds to a situation in which the 
SI of the material goes to one for every values of to, obviously an unrealistic one from the 
point of view of real materials. 

Concerning the other range: y a > 1.3, it seems to indicate that for greater and greater 
values of y a the Casimir energy gets more and more negative. This is a very delicate point 
and it deserves deeper examination. 

It is true that, in doing the calculations, we made no assumptions on the values of y a , 
but, as one can easily realize, the asymptotic expansion we did is non uniform with respect 
to the parameter y a in the range y a £ [0, oo]. This can be traced by observing that the two 
limits ya—T-oo and v — > oo do not commute, so, once we did the expansion for v->cowe 
are no more allowed to take y a as large as we want. Incidentally we note that this prevents 
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us to take the limit y a — > oo in Eq. (|2ip that would give the Casimir energy for a perfect 
conducting sphere because lim^^oo Z = 0. 

To recover this limit, instead, we had to make the limit y a — > oo first and then to 
proceed with the asymptotic expansion with respect to v. Doing in this way we find 



in agreement with |22(, 136]. We note that even in this case, the correction due to y a allows 
for a negative Casimir energy for y a < 2.3 (but, perhaps, this must be considered a too 
small value with respect to the assumption y a — > oo). 

In any case, even though it would be very difficult, maybe impossible, for real materials 
to fulfill the conditions to have negative Casimir energy, these results can be very interesting 
from the point of view of the application of the MIT bag model to the confinement of 
quarks into the hadrons. Indeed, to recover the zero-point contribution of the (static) 
MIT bag energy, E{a) ~ — obtained from the phenomenological bag model fits for the 
spectrum of hadronic particles [20j . we need simply y a ~ 1.823 (where we used Eq. (I22p 
and considered the fact that the gluons are eight) . Our result simply tells us that probably 
the hadron surface is not an ideal conductor with respect to the color. 

However a fully satisfactory application to QCD has to face with stronger problems as 
for example the self interaction of the gluons and the non abelian nature of the gauge group 
with the consequent change of the propagator of the theory [371 EH] that could eventually 
modify the result. In this respect this is only a first step. 

It would be very interesting, in our opinion, to find a uniform asymptotic expansion 
for y a G [0, oo] so to have the possibility of computing the two limits y a — > and y a — > oo 
on the same formula and to extend the range of validity of (I22p to y a — > oo too. Also 
interesting, from the physical point of view, would be to extend the same approach to 
more general Surface Impedance functional forms. 
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E, 



'cas\y a — too 



- 0.084 + 0.0081 log (fia) + 



0.0040 

s + 1/2 



+ 




(23) 
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Appendix A 



Mv) 
Ai(y) 

Mv) 
My) 



Mv) 



ViFTT y 2 + i 

y {y [y («y (ViFTT - y) - 21) + 20^+1] - s} + 4^/y 2 + 1 



8(y 2 + l) 5/2 

16y 8 + 56y 6 + 73y 4 + 10y 2 + 2 y (l6y 8 + 64y 6 + 97y 4 + 38y 2 + 4) 



8(y 2 + l) 7/2 



8(y 2 + l) z 



2y^^+l 



-80y - 512y 13 - 2720y n - 5808y 9 - 6193y 7 



128(y 2 + l) 11/2 

vV + 1 (512y 12 + 2464y 10 + 4640y 8 + 4176y 6 - 160y 4 + 288y 2 + 16) 



1520y 5 - 112y d + 
+ 



4y 2 (y2 + 1)5/2 



(-4y 7 



10y 5 



6y 3 + Vy Y +l (4y 6 + 8y 4 + y 2 - l)) 



3<5 4 



8y 4 Vy J Tl 



| 5l [{Uy s + 224y 6 + 272/ + 91y 2 + 2) y 2 ] 



Si 



16y 2 (y 2 + 1) 7/2 
-8^y 2 + 1 (8y 6 + 24y 4 + 23y 2 + 3) y 3 - 2 



+ 



(1024y 16 + 7296y 14 + 22080y 12 )y 2 



16y 2 (y 2 + 1) 7/2 128(y 2 + l) 15/2 

(36336y 10 + 33763y 8 + 8007y 6 - 1512y 4 + 3780y 2 + 72)y 2 8 



128(y 2 + l) 15/2 128(y 2 + l) 15/2 

4 (256y 16 + 1696y 14 + 4704y 12 + 6928y 10 + 5461y 8 - 42y 6 + 208y 4 + 356y 2 + 4) y 

128(y 2 + l) 7 
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Bi(y) 
B2{y) 



1 - U 2 5y 3 



2(y 2 + l) 2 8(y 2 + l) 5/2 

y (-9y 5 + 5y 3 + yftp + I (-9y 4 + 10y 2 + 4) + 12y) - 2 j 



8(y 2 + l) 9/2 2y2y^i 
l s2 ( 2y 2y i + y 2 + l\ : -25y 6 + 70y 4 - 24y 2 + 1 



-401y 7 + 928t/ 5 + 112y 3 - 112y 
128(y 2 + l) 11/2 

„y ,4/l 8y 2 + 9 \ 2 / y(y 2 -3) 8j/ 6 -21^-6y 2 -2 \ 

5980y 6 - 1363y 8 - 4292y 4 + 512y 2 - 8 1330y 7 - 341y 9 - 376y 5 - 316y 3 + 36y 
128 (y 2 + 1) 13/2 32 (y 2 + l) 7 

Appendix B 

In the following all the relevant integrations can be obtained as limiting cases of the fol- 
lowing formula 



dy 

ma/v 



(yvy 2 



y b 1 /am\b-2c+i ( 1 + 5 

' ) B [ 1 — s,c + s 



(1 + y 2 ) c 2m 2s \ v J \ 2 



1 + 6 1-6 v 2 
2-F1 c, c + s — , +c H — 



a?m 2 



with B(x,y) = ^(x+^y) ^ e beta function and 2-P1 (a, 6, c; the hyper geometric function 
[33]. Indeed one easly realizes that quite often it is possible to put m = from the very 
beginning and in this case one can use directly 

fy _vy2s yb (g)" 5 r(^)r( c + g -^) 

V \ a ) (1 + y 2 Y 2T(c) 
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In this way we find for the TE-modes: 



Co = E 

u=3/2 

oo 

= £ 

u=3/2 



i/sin(irs) f 00 / + y 2 - 1) 



+ 



2(y 2 + l) / V a 



yu 



-2s 



dy 



^sin(Trs) (f ) 2S (vr 3 / 2 csc(tt S ) - uY (s - \) T(-s)) 



TT 



4^ 



4 I \ ' 2 J 7r 3 / 2 y 



(24) 



In treating C\ we have to introduce a spurious exponent a for the y so to avoid the 
divergencies when summing over v, thus we define 



C\ = lim V 



^sin(7rs) f°° -v 2 y a 



a-^l * — ' TT 
i/=3/2 



o ^V + V 2 



a ) 



lim f; -ia 2 V 2 -« sin(7r S )j,r 2s " 1 sec (±7r(a - 2s)) 

i/=3/2 ^ ' 

lim - V( (a - 2, ^ sinM^" 2 *" 1 sec Qvr(« - 2s) 



(25) 



And for C 2 : 

sin(7rs) 



i^=3/2 



m 



A („\ , Mil , ^M. M ^ M Mvl 

Ao(y) H h - o + .q H -3— 



shifts) 9 „ 
— - — -a zs 



TT 



ft" {*) , , , /TOO 



+ 



+ 



+ 



TT 



E ^ 2s+1 

i/=3/2 

c (2, - 1, 1) / T£ ( S ) + c (2,, I) /r(*) + c (2. + 1, |) g E (s)+ 



( ( 2s + 2,1 ) / 3 ™(s) + C ( 2s + 3,^ ) fJ E (s 
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with 

fT E {s) 
fI E (s) 

fF(s) 



fF(s) 



7T 



7T 



sin(-7rs) 



a 2s ( (2s, |) (I27r 3 / 2 ssec(7rs) + (5(3 - 2s)s + 7)r(l - s)T (s - ±)) 



240F 



sa 



2s 



?r((13 - 5s)s + 16) 



32sin(7rs) 

(l2y 2 (s _ 1) + 40s 4 - 28s 3 - 58s 2 + s + 9) T (s - \) V(s - 1) 



a 2s {- 



1 



384 cos(vrs) 



480F 

7ts [96y 2 + (2s + l)(4s(15s - 7) - 169)] + 



(;o L (2s + l)(-7560y 2 (2s - 3) + s(s(s(4s(1105s - 7184) + 18431) + 

69542) - 62271) + 1764) T(l - s)T (s - 3/2) } 

2 r vr (768y 2 (2s + 3) + s(s(s((621 - 113s)s + 1691) - 2253) - 4362) - 384) 
° I 3072 sin(Trs) 

y 2 (s(s(77 - 20(s - l)s) + 4) - 9)r (-s - \) T(s - 1) 



+ 



480? 



+ 



60480^ 



s[s (s (8s (s (2210s 2 - 5968s - 10351) + 37085) + 3385) - 379346) + 
129783] + 17010)r(l/2 - s)T(s - 2)]} 



For the TM-modes we get: 
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D = 
D 1 = 



Co 

sin(7rs) 



i/=3/2 



-2s 



sm(7rs) ^ a u 



2s,,l-2s 



y 2 + i ^/yi+i 



r(i- s )r(s) 



7T — 2 

u=3/2 



sinl its 



TT 



a 2 % (2s - 1, |) (sin(7rs)r (| - s) T(s) + tt 3 / 2 ) 



2vr 3 / 2 



£>9 = 



sin(7rs) 

TT 

sin(7rs) 



V 



^ !/=3/2 • /ma ^ 



yu 
a 



— m 



Bi(y) , B 2 (y) , B 3 (y) , B A {y) 



+ 



+ 



+ 



IT 



i/=3/2 



2s+l 



l -fI M (s) + l/™(s) + l/ 3 ™(s) + 1/ 4 ™( S ) 



sin(^s) q2s 



7T 



C(2s + 3,-)/™(s) 



with 

/2™00 = 



1 



1 



5r(2 - s)r (s + 1/2) 



C0S(7TS) 

2 (3 - 12y a 2 + 2s (20s 2 + 6s + l)) V (-s - 1/2) T(s + 1) + 



96^/tt 

3tt 3/2 (13 - 5s)s 
1 



2 1 



sin(7rs) 



315vr 3/2 s (23 - 96y a 2 + 2s (60s 2 + 2s + 9)) 



COS 7TS 



+ 



120960V7T 

3 (7560y a 2 + (2s + l)(s(131 + s(1105s - 4974)) - 672)) T(l - s)T (s + 1/2) 
sin(7rs)a 2s C (2s + 3, §) r 315vr 3 / 2 s 



967680vr 3 / 2 I sin(7rs, 

s(s + l)(s(s(113s - 734) + 3) - 246) 



[768y a 2 (2s + l)+ 



4r i^-s - -J T(s + 1) ^15120y a 4 (s - 3) + 2520y a 2 (2s + 3) (20s 3 - 5s + 3) 
(s + l)(2s + l)(2s + 3)(8s(s(2s(1105s - 774) + 1363) - 177) + 945) } 
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